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I.  Introduction 


Let  w.  be  a  standard  n-dimensional  Brownian  motion,  and  let  x. 
be  an  n-dimensional  diffusion  which  is  the  solution  to  the  following 

stochastic  differential  equation: 

dx{  =  f(xt)dt  +  dw{  (1.1) 


where  fi  e  C^bCR11),  i=l,...,n.  We  are  interested  in  computing  the 
asymptotic  behavior  of 


P(ll(j)-xll<e) 

P(llwlke) 


=  J(<M); 


(1.2) 


as  £  -»  0,  where  <{>  is  a  deterministic  n  dimensional  continuous 
function  on  [0,T]  and,  for  any  \|/  e  C([0,T]  ->  Rn), 


ll\j/||  =  max  l\)/(t)l  (1.3) 

tetO.T] 

and  I  I  denotes  the  Euclidean  norm  in  Rn. 

This  problem  was  investigated  by  physicists  in  the  context  of 
statistical  mechanics  and  quantum  theory  cf.  [1],  [2].  A  rigorous 
mathematical  treatment  was  initiated  by  Stratonovich  and  carried 
out  by  Ikeda-Watanabe,  Takahashi-Watanabe  and  Fujita-Kotani,  in 
various  degrees  of  generality,  cf.  [3],  [4),  [5],  [6].  In  particular,  the 

two  last  references  treat  the  case  where  (1.1)  is  a  general  s.d.e  (i.e., 
with  state-dependent  diffusion  coefficients),  and  the  diffusion 
evolves  on  a  manifold. 

The  analysis  above  was  restricted  to  the  case  of  ((>  e  C2  (C°°  in 
[5],  [6];  however  it  seems  that  their  technique  can  be  pushed  through 
up  to  cover  C2).  In  that  case,  it  was  shown  that 


J(<)>)exp((-||0||— 1  )K(e))  <  J(<|>,e)  <  exp((4||+l)K(e)  +  |ffi|e)J(<|>)  (1.4) 


where  K(e)  ->  0  and 


3 


T  T 

J(<W  =  Urn ^  J(<J>,e)  =  exp  -[i  J(i$(s)-f(<Ks)l2)ds  +  i  JVf(0(s))ds]  (1.5) 

0  0 

In  the  context  of  the  estimation  of  trajectories  of  diffusions, 
there  was  a  need  to  evaluate  (1.2)  for  certain  <[>  which  are  not 
necessarily  C2;  for  a  specific  class  of  random  <()  (which  correspond, 
roughly,  to  <|>(t)  =  I1  vsds,  where  v.  is  a  Brownian  motion  which  is 
independent  of  w.),  it  was  shown  by  probabilistic  methods  in  [7]  that 
still  J(<j>)  =  lim  J(c>,e),  a.s.  Pv. 

Our  goal  in  this  paper  is  to  evaluate  lim  J(0,e)  for  0  which  are 
not  in  C2[0,T].  That  will  allow,  in  the  estimation  problem  considered 
in  [7],  to  include  feedback  in  the  observation  model.  The  main  result 
is  collected  in  the  theorem  below. 

Theorem  1 .  For  0e  C1+a,  a>0  deterministic,  lim  J(<|>,e)  =  J(<J>) 
where  J(0)  is  defined  by  (1.5). 

We  remark  that,  in  the  case  of  a  diffusion  evolving  on  a 
manifold  (or,  more  specifically,  in  the  case  of  state-dependent 
diffusion  coefficients),  the  functional  J(<J>)  involves  an  additional  term, 
related  to  the  scalar  curvature;  however,  the  result  J(<j>,e)  -»  J(<[>)  still 
holds,  c.f.  the  remark  in  the  end  of  section  3. 

We  note  that  Takahashi  [13,  remark  1,  page  379]  has  claimed  a 
stronger  version  of  theorem  1  and  it's  converse.  However,  no  proof 
is  given,  nor  has  one  been  published  since.  We  did  not  succeed  to 
prove  the  theorem  in  the  stronger  form  appearing  in  [13]. 

We  conclude  this  introduction  by  a  "cheap"  proof  of  our  results 
for  0  e  C1+a>  a  >  1/2,  of  a  converse  result  when  <|>  e  Aa2’°°  where 
Aa2’°°  denotes  the  fractional  Sobolev  space  (cf.  [8]),  a  <  1/2,  and  by 
some  notation  conventions.  Section  2  includes  a  description  of  the 
problem  in  terms  of  a  PDE  approximation  problem,  and  section  3 
includes  the  proof  of  our  main  theorem. 

Let  (j)  e  C1+a,  l>a  >  1/2,  and  let  0(8)  denote  the  mollification  of  <|>  by  a 
8-mollifier.  By  extending  appropriately  0(t)  for  t<0,  let  (j>(0)  =  <J>(5)(0) . 
Then  (c.f.  [8]),  110-0(8)11  <_c51+<*,  Il0(5)|l  <  c5<*-l,  and 
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P(llx-<|>||  <e)  <P(llx-<t>®||<£+c81+(X)  <  P(llwlke+5)J(<j)^)exp(K(e)(||^  <5)j|+l)+ce50t-1) 


P(ilwll<y)  =  K(y,T)  exp - —  T,  (1.7) 

(Y)2 

where  ^1  is  the  first  eigenvalue  of  the  Dirichlet  problem  in  the 
unit  ball  (c.f.[4]  and  also  below,  and  K  (y,T)Y::40K.  Therefore, 


P(IMII<£)<  J(^))exp(K(e)(nV5)ll  +  l))exPa  T(~r- 
P(llwli<e)  1  e2 


1 

(e+c81+ot) 


-)  +  ce8a  !) 


(1.8) 


By  choosing  8  =  eY,  one  gets  that  in  order  to  demonstrate  theorem  1 
we  need 


( - )e-^°  0=>y>— ^  (1.9a) 

e2  (e+ce*1+a))2  1+a 


el+7(a-l) 


e — $  i 

->  0  =>  y< - 

1-a 


(1.9b) 


and  therefore, 

3  1 

- <y<- —  , 

1+a  1-a 

and  a  solution  for  y  exists  if  a  >1/2  .  A  similar  argument  holds  also 
for  the  lower  bound,  and  the  "cheap"  proof  is  completed.  Note  also 

that  a  weak  version  of  a  converse  to  the  theorem  holds  for  <j>gAa2’°°, 
a<l/2  but  (jisAa-a'2*00-  all  a'  >  0,  where  Aa2’°°  denotes  the  fractional 
(p=2)  Sobolev  space,  c.f.  [8]:  indeed,  let  <|>  (5),  denote  the 
mollification  of  <|>  by  a  S-mollifier.  Again,  (c.f.  [8]), 


Jo  Tl$(S)|2ds  >  c^2^*1)),  ll<|)(5)ll  <  Sa-2_a';  plugging  into  (1.7),  one  has  that 


P(llx~<j)lk£)  <  C  expt-cS2^1)  +  X  T(— )  +  e5a'2-a) 
P(llwlke)  1  e2 


(1.10) 


To  show  that  the  ratio  of  probabilities  in  theorem  1  converges  to  zero 
as  e  0,  we  need  to  show  that  the  R.H.S.  of  (1.10)  0,  for  8  =  eY. 

But,  similarly  as  above,  one  gets  the  pair  of  conditions: 

1  l 

-  <  y  <  — 

1-a  a 

which  possess  a  solution  for  a  <  1/2. 

Our  goal  will  be  therefore  to  "close  the  gap"  left  by  the  cheap 
proof;  we  do  that  by  reducing  the  problem  to  the  case  of  f  =  0  (no 
drift),  following  [4],  and  then  transforming  the  problem  to  a  PDE 
one.  This  will  allow  us  to  get  much  tighter  bounds  on  the  distance 
between  the  "regularized"  solution  (with  <t>(5))  and  the  solution  to  the 
original  problem,  and  that  will  yield  the  sharp  estimates  announced 
in  the  theorem  above. 

Notations 

Throughout,  Q.  denotes  the  unit  ball  in  Rn,  and  e£2  denotes  the 
ball  in  Rn  with  radius  e.  II  Ilk  denotes  the  k-th,  p=2  Sobolev 

norm  in  Q,  i.e. 

ii(t>llk  =  (  ^  J(Da<j))2dx)1/2 
|a|<k 


where  the  domain  of  the  integration  (£2,  eQ.)  will  be  clear  from  the 
function  involved. 

v*  denotes  the  transpose  of  a  vector  v. 

u*v  denotes  the  composition  of  u  and  v,  c.f.  section  3. 

II  II  denotes  the  sup-norm,  and  I  I  denotes  the  Euclidean  norm 
in  Rn. 
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II.  An  associated  PDE  formulation. 


In  this  section,  we  reformulate  (1.2)  in  terms  of  an  associated 
PDE.  A  similar  approach  can  be  found  also  in  [5]. 

We  start  by  noting,  following  [4],  that,  for  <j)eC1[0,T],  xt~<j)(t) 
satisfies: 

d(xt-(>(t))  =  $(t)dt  +  f(xt-<|>(t)  +  <Kt))dt  +dwt  (2.1) 

By  Girsanov's  transformation,  one  has: 

T  T 

P(llx~<i)|l<e)  =  E(exp(  f(f*  (w  +<|>(t))  )dw M  If  «J)(t)+wr)-«5(t)l2dt)|llwli<e) 

P(llwlke)  o  o 

(2.2) 


Note  that 

T  T  T  T 

Jf  (wt+<j)(t))dwt  =  Jf  ((<j)(t))dwt  +  JwtVf*((|)(t))dwt  + J  0(w2)dwt, 
ooo  o 

where  Vf  denotes  here  the  matrix  of  partial  derivatives  of  f  and 
also 

T  T 

r  f  X"'  r)f®  j 

j  f*(<>(t))dwt  =  wTf*(4)(T)) 2-i  wt  ^dt  (2-3) 

o  o  1,J  -> 


and,  by  Ito's  lemma 


f  *  IwJ2-T  J  (IwtI2-T) 

J  w*  Vf(4)(t))dwt  =  V*f(<>(t))(— T— )  -  J  -----  ( V  V  .f(0(t))  4>(t)dt 


+  J  X(— }  (<Kt))wtW  =  (lljll  +  H<t>H)0(e2)  +  j  Jv  •  f(<|>(t))dt 

0  i*j  3Xj  0 


+ 


J  X  (— 1)(<t>(t))w1tdwJt 


(2.4) 


o  *j  9xj 


where  V*f  denotes  here  the  divergent  of  f;  combining  (2.2),  (2.3) 
and  (2.4),  one  has: 

P(-'-X--)  =exp(-y  [ lf«i>(t))  -  J(t)|2dt  -  J  fv  •  f((Kt))dt) . 

P(llwlke)  1  i 


E(exp<0(e2)(|i(f>!|  +  ||(|)||)  + 


J  2L  (4)(0)w1dw 

o  *idxi 


T  T 

+  J  0(lwl2)dwt  -  J  $*(t)dwt)l  llwll  <  e) 
o  o 


(2.5) 


By  lemmas  of  [4,  pg.  451]  (see  also  [6]),  (which  are  the  main  part  of 
the  proof  in  [4]), 

T 

E(e?p  c  j*0(lwl2)dwtl  llwlks)  £:401  Vc  (2.6a) 

o 

T 

E(expcJ  ky^w'dwj  |  HwII<£)£;:401  Vc  (2.6b) 

o 

and  therefore,  to  compute  (2.5)  we  need  only  compute 
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T 

E(exp  -J<j>*(t)dwtl  llwlke) . 
o 

and  show  that  it  converges  to  1  as  e  ->  0.  Let  us  define 

T 

A  =  E(exp  -  J  }*(t)dw  J  llwlke)P(llwike). 
o 

Then,  by  Girsanov's  theorem 

T 

A  =  exp-i(Jl$(t)l  dt)P(ilw-(J>||  <  e) 
o 

Let  u(z,t,x,s)  be  the  fundamental  solution  of 

1  12 
Ut  =  -  Au  +  <5(t)Vu  +  —\  <j>(t)  I  u 


u(z,t,x,s)  |  =  0  (2.7) 

izl=£  V  ' 

i..e.  the  solution  of  (2.7)  such  that,  for  each  continuous  f(x), 

lim  u(z,t,x,s)f(x)  dx  =  f(z) 

t: — >s  ^  ^ 

eQ 

Such  a  solution  exists  and  is  unique  by  the  maximum  principle 
(cf.  [9,  ch.  1-2]).  Then 

A  =  Ju(z,T,0,0)dz 

e£2 


(2.8) 


Our  goal  will  therefore  be  to  compute  bounds  on  the 
fundamental  solution  of  (2.7).  It  turns  out  that  one  can  find 
explicitly  the  solution  to  a  related  equation  (eq.  2.3a),  and  then  by 
perturbation  techniques  relate  the  two.  Towards  this  end,  let  <|)(8) 
be  a  8  mollification  of  <(>  (for  example,  with  a  Bessel  potential,  or 
otherwise,  cf.  [8]),  and  let  u(§)(z,t,x,s)  be  the  fundamental  solution 
of: 


u®  =  J  Au(S)  +  $(t)Vu(5)  +  ($(5)(t)  -  $(t))Vu(5)(t)u(5> 


(2.9a) 


U(8)(z,t,x,s)|  =0 

|zl=e 


(2.9b) 


In  the  sequel,  let  j(8)(t)  =  <|>(8)(t)  -  <|>(t).  We  will  assume  throughout, 
without  mentioning  it,  that  ||j(8)(t)||  <  1. 

Our  line  of  attack  will  be  as  follow:  we  first  show  below  that 


f  — (z,t,0,s)  £-*0 j  fort_s>£  >  0  uniformly  in  8>0,  i.e.  that  if  in  (2.8)  one 

J  p(llwlke)  0 


eQ 


substitutes  u(8)  instead  of  u  one  has  the  required  convergence 
(lemma  2.2).  We  then  show  in  section  3  that 


jdz(u®(z,t,0,s)  -  u(z,t,0,s)) 
eQ _ 


P(llwlke) 


e— >0 
->  0, 
8(e)-»0 


where  8(e)  -»  0  in  an  appropriate  way,  thus  establishing  the  required 
convergence.  To  demonstrate  this  last  convergence,  note  that  the 
solution  to  (2.7)  can  be  represented  by  the  classical  parametrix 
method  in  terms  of  an  infinite  series  involving  the  solution  of  (2.9) 
(theorem  3.1).  Estimates  on  u8(z,t,x,s)  which  we  prepare  in  the 
remaider  of  this  section  are  crucial  in  obtaining  the  required 
convergence. 

We  use  the  following  classical  result: 


Lemma  2.1 

u(5)(z,t,x,s)  exists  and  is  unique.  Moreover,  there  exists  a  c 
independent  of  e,  5,  such  that 


(8)  C 

iu  (z,t,x,s)l  < — ~2  exp 

(t-s) 


(  2 
(z-x) 

C(t-S)  j 


(2.10a) 


cS’s  c 

IVu  ;(z,t,x,s)l  < - ^7-  exp 


f  2) 
(z-x)25 


(t-s) 


n+1 

2 


(  c(t-s) 


In  particular, 


(2.10b) 


IVu  ®(z,t,x,s)l  < 


(t-sflz-Xr2^ 


V  1/2  <  p  <  1 


(2.10c) 


IVu^(z.t.x.s)!  < - - 


(t-s)^lz-xln+1'2^ 


(2.10d) 


Proof.  The  estimates  (2.10a)  and  (2.10b)  are  the  well  know  Arronson 
estimates.  For  an  easy  derivation  of  them,  we  refer  to  [12]  and 
references  there.  (Note  that  in  [12],  only  (2.10a)  is  proved,  however 
(2.10b)  follows  easily  by  differentiating  throughout  in  the  proof). 
(2.10c),  (2.1  Od),  which  are  the  only  estimates  we  will  need,  follow 
easily  from  (2.10a),  (2.10b).  A  different,  more  cumbersome  proof  of 
(2.10c),  (2.1  Od)  via  the  parametrix  method  appears  in  [9,  ch.  1, 
section  4-5].  Finally,  uniqueness  follows  from  the  maximum 
principle. 

The  usefulness  of  equation  (2.9)  lies  in  the  fact  that  it's 
solution  is  easily  represented;  to  do  that  we  need  some  auxilliary 
results,  which  are  regrouped  in  (a)-(c)  below;  (2.13)  is  the 
representation  of  the  solution  we  will  use  in  the  sequel. 

Lemma  2.2 

Let  (ym(x),  Xm)  denote  the  normalized  (w.r.t.  L2(£2)) 
eigenfunctions  and  eigenvalues  of  the  Dirichlet  problem  in  the  unit 
ball  in  Rn,  i.e. 


1  2 


Aym(x)  = -^mym(x),  x  s  n  (2.iia) 

TmW  1^1  =  0  (2.11b) 


Then: 

(a)  There  exists  a  unique  eigenvector  associated  with  the 
minimal  eigenvalue  Xo,  and  Xq  >  0. 

(b)  The  set  {km}  is  discrete,  and,  if  N(A,)  denotes  the  number  of 
eigenvalues  s.t.  Xm  <  X  (including  multiplicity),  then 

N(X)  =  KXn'2  +  oQi^2),  as  X  ->  ~  (2. 12) 


(c)  Ym(x)  e  L2(Q);  Moreover,  Vk,  ||Ym(x)l|k  <  00  and  |Ym(x)|  < 

c(?im)n/2: 

Finally,  Ym(x)  spans  1_2(£2). 

(d)  The  following  limit  exists  (pointwise,  uniformly  in  (z,x) 
for  t-s  >  xoe2,  for  all  e  and  in  L2(eOx[0,T]))  and  is  the  fundamental 
solution  of  (2.9): 

(ft  ('S')  A  x1  1  A,  (t-s)  Z  x 

u  (Z,t,x,s)  =  lim  u:  (z,t,x,s)  =  lim  2*  exp  -( - -)Ym(~)Ym(~) 

j— »<x>  j— m=0  £  £^  £  £ 


L 

exp(-($  %)(z)  -  $  ®(s)x))  exp  (i  J  l<J><8)(x)  -  <J>(x)/2dx  (2.13) 


Proof:  For  (a),  c.f.  [10].  (b)  is  theorem  (14.6)  of  [11].  That  Ym(x)  e 
L2(^)  follows  from  theorem  (16.5)  of  [11].  To  see  that  ||Ym(x)||  < 
note  that 


"Avm(x)ii0  =  amuiA‘Vm(x)ii0  =  ajt 


(2.14) 


Therefore,  by  the  Sobolev  lemma  (c.f.,  e.g.,  theorem  (3.8)  of  [11] 


Moreover,  since  WQln/2]+1  C*  C(£l),  one  has  also 


([~]+ 1/2) 

ly  (x)l  <  kA, 

‘mN  7  m 


(2.16) 


Finally,  we  show  (d).  Note  that 


£ 

m=ll 


1  \,(t-s) 

—  exp  - 


Ym(-)Ym(~)  exp  -(J[5)z4(s5)x)exp(y  f Ic|j(6)(x)  —  <p(T:)[2dx) ! 

p  p  ^  •/ 


£4 


<  /  .  —  exp 
m=i  e  e 


,  (r 11 


„ 

2  nj 


.(f +1)  -(M-)  ^ 

j  exp  e2  '  <  00 


(2.17) 


where  we  have  used  (2.12);  note  that  the  convergence  is  uniform  for 
t-s>e2xo,  is  independent  of  e  for  t-s>e2xo  and  also  that  it  holds  even 
after  scalling  by  exp  -(Ao(t-s)/e2).  The  convergence  in  L2(eQx[0,T]) 
is  very  similar,  and  will  not  be  used  in  the  sequel. 

It  is  easy  to  check  ,  similarly,  that  the  convergence  holds  also 
for  the  derivatives  of  Uj(5)(z,t,x,s)  (w.r.t.  t  (once)  and  w.r.t.  z  twice), 
and  that  lim  Uj(8)(z,t,x,s)  satisfies  (2.9).  It  remains  to  check 
therefore  that  it  is  indeed  a  fundamental  solution. 

Let  f(x)  be  a  Co°°  (on  eQ)  function  (and  in  particular, 

Yi(— )£  with  2jl  <  oo).  Let 
i=l  £  i=l 


@(t,s,z)  =  Ju^(z,t,x,s)f(x)dz- f(z) 


(2.18) 
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ll©(t,s,z)lP  <  Z.  f?[._  x 

0,£fl  1  *  2 

1=1  e 


X.(t-s) 

exp  — L_  - 1] 


Let  k0  be  such  that 


i=k„ 


and  xo  such  that 


\\ 

exp  ( - )  >  l  -  v 

o  * 


one  has  then 


H0|f  <  2y 
"o.e.n  1 


for  t-s<xo,  and  since  y  is  arbitrary,  we  have 


Men) 
©  ->  0 


Similarly, 


k  ^2  -X-Ct-s) 

HA 0(t,s,z)llo ^^2,  gf(k)  [exp  — ^ - 1] 

i=l  £ 


where 


Akg(x)  =  ^  y.  (-)gj(k)  and  ^  gf(k)  < 


i=i  e 


(2.19) 


As  above,  one  has  then  that  ||Ak  0(t,s,x)||o,ef2  -»  0,  which  implies  by 
the  Sobolev  lemma  that  @(t,s,z)  0  pointwise.  Therefore,  one  has 

that,  in  the  sense  of  distributions  in  D'(e£2,  lim  Uj(8)(z,t,x,s)  is  equal 
to  the  (unique,  by  [9,  ch.  2])  fundamental  solution  of  (2.9).  Since,  as 
is  easily  checked,  for  (t-s)  >  0  both  this  limit  and  the  fundamental 
solution  are  continuous  in  z,x,  they  are  equal  everywhere,  which 
concludes  the  demonstration  of  the  theorem. 

We  establish  below  some  estimates  which  will  turn  out  to  be 
useful  in  the  perturbation  analysis  of  section  3; 

Lemma  2.3 


U(5)(z,t,x,s)  =  —  exp 0- [exp(-(c|>®(t)z-<j>(8)(s)x))]7 ’  (~)  Y0(— )  +  A(z,t,x,s)] 

e11  e2  0  e  e 

(2.20a) 


where 


k  exp  —ft'1 -  S--  if  (t-s)  >  xne2,  AX  —  X  -X 
c2  u  10 


IA(z,t,x,s)l  < 


ken 

-  if  (t-s)  <  x  e2,  1/2  <  \i  <  1 

.  ,n-2p  o  ^ 

(t-s)^lz-xl 


(2.20b) 


and  k  is  independent  of  e,8. 

Similarly, 

1  ^  ^  Z  X 

Vu(8)(z,t,x,s)  = - -exp  — 4 —  [exp  -($  (5)(t)z  -  $  (5)(s)x(Vy  (-)y  (— )  - 

en+1  e2  0  e  0  e 

-  $  (5)(s)e  y  (~)y  (-))  +  B(z,t,s,x)] 
u  £  0  £ 


where 


(2.21a) 


1  6 


kexp 


-AA,(t-s) 


IB(z,t,x,s)l  <  • 


ke 


n+1 


(t-s)^lz-x|n+1'2H 


if  (t-s)  >  T0e2 


if  (t-s)  <  tQe2, 1/2  <  (i.  <  1 


(2.21b) 

Proof.  The  upper  bound  in  (2.20b)  and  (2.21b)  follows  immediately 
from  the  representation  (2.13)  and  the  method  of  proof  of  lemma 
(2.2).  The  short  time  estimates  (the  lower  line  of  (2.20b)  and 
(2.21b)),  follow  directly  from  the  derivation  of  [9,  ch.  1,  section 
3,4],  or  from  [12]. 


satisfy 

t-S  >  £2T 

0 

0  <  t  -  s  <  Toe2  1  >  (i.  >  0,  n+1  >  n+Yj  >  0 

(2.22) 


Lemma  2.4 

Let  Cj(z,t,x,s),  i=1,2 

!,  -AX(t-s) 

k.  exp - 

1  £2 
k— £5— 

1  /  ,  ,n^i 

(t-s)  Iz-xl 

Then 


iJJ  C.(z,t,x,s)dxdsl  <  k 


se£2 


,  (n+Pryi+2(l-Ri))A(n+2) 

k.£ 

1 


(2.23) 


(with  similar  bound  when  the  integration  w.r.t.  x,s  is  replaced  by  an 
integration  w.r.t.  z,t),  and 


|C.*C(z,t,x,s)|=  ij  J  C^ZjtjX'jSOCCx'jS'jXjSjdx'ds'l  < 


t  e£2 


kk  k  e<n+<° A Pi-7i+2( 1  -^i)AOj-Yj  +  2(Hij ))  exp  -A^t-s)...  t-s  >  e2x 


i  j 


kk.k.  e 


Pi+Pj+n-VYj-2^-2^+2 


,  t  -  S  <  e2x 

7  n 


t-s  ovCtij+^-l)  z-x  ov^+yj+n) 

(  )  ] - 1 

e  e 


(2.24) 


Proof 


First,  note  that 


t 

J  J  IC(z,t,x,s)ldxds  <  k.  J  J  —  ^ ^  dsdx  +  k.  J  J 

s  z£l  s  e£2  t-x  e2  sQ 

0 


n+Pi 

e _ : 


dsdx 


Hi  n+y- 

(t-s)  Iz-xl 


<  (kk.e1*2  +  kk.e 


n+Pi+ZCl-jLtjJ-Yi 


from  which  (2.23)  follows. 

Considering  (2.24),  let  first  (t-s)  <  e2xo;  one  has  then 

t 

|I1(t— s)|  =  JJ  |Ci(z,t,x',s')NCj(x',s,,x,s)ldx’ds  < 

$eQ 


L. 

£kikJ 


ds’ 


J 


^+^+11-7^ 


-dx' 


-  ,  Pi,  ,  Mj  J  ,Z  n+Yi  ,  X  n+Tj 
s  (t-S)  (S  -S)  Q  \—  -  X  I  lx - I 


(2.25) 


We  recall  the  following  (c.f.,  e.g.,  [9,  pg.  14]): 


1  8 


J 


dx' 


|a-x1a|x’-b|p 


< 


k  la-bln"a  P 
k 


if  n  <  a  +  p 
if  n  >  a  +  P 


Applying  (2.26),  one  has  from  (2.25) 


(2.26) 


I  I1(t-s)l  <kk.k. 


I 


(\+Pj+2n 

e 


Iz-x 


„n+Yi+Tj 


( if  Yi+Yj+n>0,  ji+pM,  t-s<e2to) 


(2.27) 

with  similar  bounds  for  the  other  cases  of  t-s<e2T:o-  We  consider 
therefore  now  t-s>£2to.  In  this  case,  we  get: 


s+e  x. 


o 


II1(t-s)l 


J  J* 


-AX(t-s') 


n+pj 


exp 


efi 


Rj  n+7j 

(s'-s)  V-xl  J 


dx'ds' 


t-£  X 

+  J  J  k.k.  exp  :AX(^  dx'ds' 

2  „  £ 

s+n  e£i 


t-e2x0 


M-i  n+y. 

(t-s’)  Iz-x’l 


Using  (2.22),  one  gets 


(2.28) 


IL(t-s)l  <  kk.k.  exp  W-V) 

*  1  J  n 


<  tt.k.  exp  (en+t°A(Pi  •  r,  -  *  WP,  -  is  *  ^  (2.29) 

1  J 


and  the  lemma  is  proved. 


III.  A  Solution  to  the  Original  PDE 


In  this  section,  we  construct,  by  a  perturbation  method,  a 
solution  to  eq.  (2.7),  based  on  u(8)(z,t,x,s). 

Let  L  denote  the  operator: 

L  v(z,t,x,s)  =  -j(5)(t)*Vv  (z,t,x,s)  +  z<[)5)(t)v(z,t,x,s)  (3.1) 

t,0 


As  before,  let  *  denote  the  composition  of  two  functions  in 
the  following  form: 


t 

f^z.t.x.s)  *  f2(z,t,x,s)  =  J  J  f1(z,t,x',s’)f2(x',s',x,s)dx<ds 

s  f.q 


Define 


Lj  §  u(8)(z,t,x,s)  =Lt  5u(5)(z.t.x.s) 


(3.2) 


(3.3a) 


and 


Lk  u(5)(z,t,x,s)  =  L  u(6)  *  Lk  1  u(5) 

t,o  t,o  t,8 


(3.3b) 


Let 


uj(z,t,x,s)  =  u(S)(z,t»x,s)  +  ^  u(8)(z,t,x',s)*L^  §u(5)(z,t,x,s)  (3.4) 

i=l 


Finally,  assume  that  e||^5)||<1  (which  is  possible  if  5  is  chosen  not 
too  small).  We  will  show  that: 

Theorem  3.1 

(a)  For  any  (t-s)  >  0,  ui(z,t,x,s)  converges  (uniformly  in 
z,x  e  eQ)  to  a  limit  u(z,t,x,s); 

(b)  u(z,t,x,s)  is  the  fundamental  solution  of  (2.7). 

(c)  Let  y(e)  =  e||(])®||e2|00,  and  let  llj®(t)ll  =  0(e%)  f or  some  X>0  . 


Then,  for  any  to>0  and  (t-s)>xo, 


(3.5) 


exp  — - En  I  u(z,t,x,s)  -  u®(z,t,x,s)l  E;40  0 

e2 

uniformly  in  z,s,  e  eQ,  and  the  rate  of  convergence  is 
controlled  by 


(5)  iffUS)  -J&lli-H  AA/c 

£||$  ||  + exp  (-? - ™ - H  1/2  <  |j.  <  (l+%)/2  (3.6) 

el/(l+)  g2 

Proof 

Part  a:  Note  that,  by  lemma  (2.3), 

-Vt-s) 

exp  - 

L  u(5)(z,t,x,s)  =  exp(-($5)(t)  z  -  $(5}(s)x)) - ^ —  [a  (t,z)y  (-)y  (-)  + 

t,8  Pn  1  0  e  u  e 


z  x  -A,  (t-s) 

+  P.  (t)Vy  (~)y  (-)]  +  exp  — & - E  (z,t,x,s)  (3.7) 

1  0  e  0  e  e2  1 


where 

ldi|  <  e||$(S)||  +  ||V5)||  ||j(8)ll  =  l(S)e4°  0  (3.8a) 


j8(t) 

IP  l  =  li]l 
1  e 


(3.8b) 
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Kk 


(6) 


(t-sflz-xln+1'2R 


t-T  <t£2 
0 


lEjl  < 


K  e*p(:^V> 


^n+1 


where 


t-s  >  x0e2 


(3.8c) 


k(5)  =  llj(5)ll  +  e2||(jj®||  e4°  0 

by  our  assumptions. 

Using  lemma  2.4,  one  obtains: 


— ) 


t,8  t,8  n 


Z  X 


Z  X 


•exp  -($  (5)(t)z-$  (S)(s)x)[d  y  (~)y  (-)  +  p  Vy  (-)y  (-)+y  ] 

ZUpUp  2  UpUp  2 


£  u  £ 


£  u  £ 


-X  (t-s) 

+  exp  — y - E2(z,t,x,s) 


(3.9a) 


where 


ld2|  <  (K'K1(5))2 

(3.9b) 

•S 

ipj  <  (K'H1(5)ii—  II 

o 

(3.9c) 

£ 


ly>|  <  Hj5(t)IIK'k(5) 


(3.9d) 
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IE2(z,t,x,s)l  < 


K'Kk 


(6)  exp 


-AA,(t-s) 


(K'Kk®)2 


(t-s)2H"1lz-xln+2(1'2fi) 


t-S>£2T 

0 


t-S<£2T 

0 


(3.9e) 


Since  2>2|i>1,  the  singularity  in  (3.9e)  is  weaker  than  that  of 
(3.8c).  By  the  same  reasoning,  one  obtains  that  there  exists  a  ko 
such  that 


u5(z,t,x,s)  =  —  exp  -(<|>®(t)z  -  <|)®(s)x)exp  ■)  [a,  y  (-)  + 

pn  P2  *o  ®  £ 


+  exp  exp  ( - 


(3.10) 


where 


la.  |  <  (K'Kl(5))k°  0 


(3.11a) 


l{3,  |  <  m®(K'K) 


.(5) . . 

,xrv  0 14 — L||^  ^th  ||m®llE::400,  at  least  as  fast  as  k® 


(3.11b) 


lv,  |  <  n®(K'K)k°  with  I  In®  1 1  £=^  0  at  least  as  fast  as  k® 
*0 


and 


(3.11c) 


(3.1  Id) 

Therefore,  by  the  same  argument  as  in  [9,ch.1,  eq.  4.7],  one  has  that, 
for  k>ko 

Lk  (u(S)(z,t,x,s) )  =  —  exp  -(V8)(t)z  -  V5)(s)x)exp  — —  -{a.  y  (-)  + 
t,5  £n  e2  k  o  e 


~  z  -A,  (t-t)  -AX/t-s') 

+  Pk  Vy  (-)  +  y  k]  +  exp  (-2- — )  exp  ( - -)  Ek 


(3.12) 


with 


lak|  <(K’Kl(5))k 
.(8). . 

Ipkl  <  (l(5))k  11^— ^-llk  (K’K)k 
e 

lyj  <  (llj5(t)IIK'Kk(5))k 


(K’K(t-s)  ^)  k 5  k  ! 

IE,  I  <  - - —  — ■  ( - )k  — 

k  T((l-|i)k+l  e  en 


(3.12a) 


(3.13b) 

(3.13c) 


(3.13d) 


By  computing  u(8)(z,t,x,s)  *  Lkt,s  (z,t,x,s),  the  Pk  term  drops  out  due 
to  the  integration 


one  has  therefore,  for  t-s>0: 
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iT  u(8)  *  Lj  u(5)(z.t.x.s)l 

Mmm  1,0 

j=k 


<  Kexp 


£ 

i-i1 


( a(e)1  +  exp 


-AX(t-s) 


r((i-p)j+i) 


e 


uniformly  in  z,x. 

(Remark:  a  similar  proof  holds  also  for  the  first  two  z 
derivatives  of  u,  for  details,  of.  e.g.  [9]). 

Part  b:  The  proof  is  identical  to  the  one  given  in  Lemma  (2.2),  due 
to  the  fact  that  u(5)*Lt,sku(8)  has  a  weaker  singularity  in  the  origin 
than  u(5);  we  ommit  the  details. 

Part  c:  By  (3.10)  (3.13)  and  the  fact  that  by  comparing  with  u<5), 
the  j3|<  term  drops  from  (3.10),  (3.12),  one  has  that,  for  t-s  >  to.  and 
e,  8  small  enough, 


lu(5)(z,t,x,s)*Lk tu(5)(z,t,x,s))l  <— 


t,S 


.  J.  (8)  ,  (5).  k 

<  —  (K(k  +1  )  exp  ■  0 


+  exp 


-\(t-s)  -  AA,(t-s)  r/j®K  yk  1  -i 

~~ —  exp  [{+ - )  • - J 


r((l-p)k+l) 


Therefore, 


en|u(z,t,x,s)  -  u(5)(z,t,x,s)lexp  — ^ [(K  (k(5)  +  l(5)))k  + 

e2  k=i 


.•(8)t 


Ad-n) 
vAX(t-sV  1  P  ' 

exp  (_7_)  (_^i — >  - 
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<  (1(5)  +  k(5)  )K  +  exp(- 


-AXT  Ki(5)  77  i(5) 


(3.15) 


By  our  assumptions, 


=  0(ex-1),  y> 0; 


Let  1/2  <  |i  <  (1+%)/2.  In  this  case,  the  R.H.S.  of  (3.15)  is  bounded  by 

-AX.X  tr 

(lw  +  kW)  K  +  (i— )  exp  ( - fl)exp  +  (- - -),  q>0  (3.26) 

e  e  e2-q 

Therefore,  (3.5)  holds  and  moreover  the  rate  of  convergence  is 
controlled  as  in  (3.6). 

Corollary  3.1.  Assume  <)>  e  C1+a,  a>0;  then  theorem  1  in  the 
introduction  holds. 

Proof.  Let 

0  <  y  and  8  =  eY ,  then  lle^5)||  < k  e  8a~1£:400 

l-a 


and 


llV5)-4ll  <  KeaY, 

so  that  the  conditions  of  theorem  (3.1)  hold.  By  (2.5)  and  (2.8),  one 
has 

J(<|>)  exp(0(e2)(||^||+||5»||)+K(e))  f  u(z,t,x,s)dz 

P(Hx-(|)||<e)_ _ gO _ 

P(llwlke)  P(llwlke)  ^ 


J(<|))exp(0(e2)(||<]>|MI4>||)  +  K(e))  J  u(8)(z,t,0,0)dz 

_ eQ _ 

P(llwlke) 


m  exp  (0(£2)(||4>||+||$||)+K(8))  J(u(z,t,0,0)  -  u(8)(z,t,0,0))dz 


P(llwlke) 

(3.17) 

Note  that 

Ju(5)(Z,t,0,0)dz  Jy0(z)dz  yo(0) 

^ - =  1 

P(llwlke)  J y0(z)dz  yQ(0) 

a. 

combining  (3.17)  with  theorem  (3.3)  yields  the  corrollary. 

Remark.  We  remark  briefly  on  the  case  of  diffusions  evolving  on 
a  manifold  (or,  more  specifically,  diffusions  with  state  dependent 
diffusion  coefficients). 

In  that  case,  [6]  have  proved  that 

T  T  T 

T .  i  .  rl  f  i$(t)  -  f(<Kt)i2dt  1  f  .......  1 

J(<|>)  =  exp-[-J - - - +-J  div fC<|»CO)dt  +  j2-JR(<Kt))dt] 

ooo 

where  R(x)  is  the  scalar  curvature  at  the  point  x  and  the  divergent  is 
taken  on  the  manifold,  we  refer  to  [6]  for  the  definitions  involved, 
and  we  point  out  where  [6]  used  the  assumption  that  <>(t)  existed  and 
how  that  can  be  avoided. 

We  recall  some  notations  from  [6]:  a  system  of  normal 
coordinates  is  defined  around  4>(t) ,  and,  in  this  system, 

aij(t,x)  =  8ij  +  (t,0)xmx1  +  0(x3)  (3.18) 


and  we  define  the  y  process: 
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dyt  =  ®ik(t.y1)dWj  +  y1(t,yt)dt  (3.19) 

k=l 

where  |y'(t,yt)|  =  0(yt)  and  the  exact  form  of  y  is  unimportant  to  our 
current  needs.  We  recall  again  that,  by  [4,  pg.  451],  if  |qs|  =  0(e2) 
under  the  conditioning  ||w||  <  e  is  an  adapted  process,  then 

T 

E(exp  c  J  qsdwsl  lwsl  <  e)  1,  Vc  .  (3.20) 

o 

Referring  now  to  the  proof  in  [6],  we  note  that  the  only  place 
one  needed  the  existence  and  boundedness  of  <j>  was  while  attempting 
to  use  theorem  2.2,  pg.  442:  using  their  notations,  one  has  to 
compute 

A  =  E(exp  c  J*Fi(t)  d\vj  llyll<£  (3.21) 


where  F(t)  =  f(t,0)  -  $(t),  and  f(t,0)  is  in  C1  w.r.t.  t.  Assuming  also 
that  $  is  in  C1,  [6]  used  the  following  estimate: 
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T  T 

J  F*(t)dwt  =  J  F*(t)a’1(t,yt)(dyt  -  y(t,yt)dt) 
o  o 


T  T 

=  J  F*(t)a1(t,yt)dyt  -  J  F  (t)a1(t,yt)y(t,yt)dt 
o  o 

(3.22) 

since  |y(t,yt)|  =  0(|yt|)  =  0(e)  under  the  conditioning,  the  contribution 
of  the  second  integral  is  negligible.  Considering  the  first  integral, 
note  that  in  the  normal  coordinates, 

cf1(t,yt)  =  l  +  0(yj)  (3.23) 


By  (3.20),  the  contribution  of  the  second  term  is  again 
negligible,  and  therefore  we  are  left  with 

T 

E(exp(c  J  F*(t)dyt)  I  llyll<e) 
o 

In  the  case  that  F*(t)  is  C1  (which  result  from  the  assumption 
4>(t)  e  C1),  integration  by  parts  yields  the  pathwise  convergence 
(under  the  conditioning  ||y||<e).  In  the  general  case,  however,  (3.21) 
reduces  to  show  that 


X 

B  =  E(exp(c  J"(|)*(t)dy  t)  II  y  II  <  e)£— >°1,  Vc 


where 


(3.24) 


dy  =  (I  +  c(y  t) )  dwt  and  c(y  t)  =  0(e2) 


The  procedure  which  led  to  our  estimates  for  the  case  c  =  0 
can  be  repeated,  where  now  the  operator  L  includes  in  it's  first  order 
term  an  additional  term  of  the  form  k  e2(<j>(5))2,  which  turns  out  to  be 
negligible.  There  is  an  even  more  direct  way  to  see  that,  based  on 
lemma  (2.1)  of  [6]  or  again  on  a  version  of  3.20  (c.f.  [6],  pg.  449];  we 
ommit  the  details  here. 
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